Universal phase shift and non-exponential decay of driven single-spin oscillations 
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We study, both theoretically and experimentally, driven Rabi oscillations of a single electron spin 
coupled to a nuclear spin bath. Due to the long correlation time of the bath, two unusual features 
are observed in the oscillations. The decay follows a power law, and the oscillations are shifted 
in phase by a universal value of ~ n/4. These properties are well understood from a theoretical 
expression that we derive here in the static limit for the nuclear bath. This improved understanding 
of the coupled electron-nuclear system is important for future experiments using the electron spin 
as a qubit. 
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A quantum bit is engineered such that its coupling to 
the disturbing environment is minimized. Understanding 
and controlling this coupling is therefore a major sub- 
ject in the field of quantum information processing. It 
is not solely the coupling strength but also the dynamics 
of the environment that governs the quantum coherence. 
In particular, the limit where these dynamics are slow 
compared to the evolution of the quantum system is in- 
teresting. The well-known Markovian Bloch equations 
that describe the dynamics of a driven system, including 
the exponential decay of the longitudinal and transverse 
magnetization then lose their validity. Such devia- 
tions from the exponential behavior have been studied 
theoretically 0, IE] and experimentally, for instance in 
superconducting qubit systems 

An electron spin confined in the solid state is affected 
predominantly by phonons via the spin-orbit interaction 
|E IE El IE HI > and by nuclear spins in the host material via 
the hyperfine interaction. At low temperature, coupling 
to the nuclear spins is the dominant decoherence source 
0, EH, El EE 0, EE EE Although this strong 

coupling leads to an apparent decoherence time T 2 * of 
the order of 20 ns when time-averaged over experimental 
runs, the decoherence time T 2 strongly depends on the 
dynamics in the nuclear spin bath. This typical nuclear 
spin dynamics is very slow, because the nuclear spins are 
only weakly coupled with each other and the bath itself is 
coupled very weakly to its dissipative environment (like 
phonons). This implies that here, the Markovian Bloch 
equations are not valid. 

Here we study the dynamics and decoherence of an 
electron spin in a quantum dot that is coherently driven 
via pulsed magnetic resonance, and is coupled to a nu- 
clear spin bath with a long correlation time. We find 
experimentally that, remarkably, the electron spin os- 
cillates coherently, even when the Rabi period is much 
longer than T 2 = 10 — 20 ns. In addition, the character- 
istics of the driven electron spin dynamics are unusual. 
The decay of the Rabi oscillations is not exponential but 
follows a power law and a universal (parameter indepen- 
dent) phase shift emerges. We compare these experi- 
mental results with a theoretical expression, derived in 



the limit of a static nuclear spin bath. 

We consider a double quantum dot with one electron 
in each dot and a static external magnetic field in the 
z-direction, resulting in a Zeeman splitting e z = g(isB z . 
The spin transitions are driven by a burst of a trans- 
verse oscillating field along the rr-direction with ampli- 
tude B ac and frequency to, which is generated by a cur- 
rent I s through a microfabricatcd wire close to the double 
dot [13]. The interaction between the electron spin and 
the nuclear bath is described by the Fermi contact hy- 
perfine interaction S ■ h, where h is the field generated 
by the nuclear spins at the position of the electron. For 
a large but finite number of nuclear spins (N ~ 10 6 for 
lateral GaAs dots) h z is Gaussian distributed (due to 
the central-limit theorem ) w ith mean ho = h z and vari 
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12| . For a sufficiently large 
3> c), we may neglect the 



external magnetic field (e 
transverse terms S± ■ hj_ of the hyperfine interaction that 
give rise to electron- nuclear-spin flip-flops (sec below). 
Furthermore, if the singlet-triplet energy splitting J is 
much smaller than both e z and g(iBB ac , we may treat 
the spin dynamics of the electrons in each dot indepen- 
dently (valid for times less than 1/ J) . 

For each dot we thus have the following spin Hamilto- 
nian (h~ = 1): 



H(t) = -(e z + h z )a z 



■ cos(tot)a x , 



(1) 



where cr, (with i = x, z) are the Pauli matrices and 
b = g/isB^c (taken to be equal in both dots), lon- 
gitudinal nuclear field h z . Here, h z is considered as 
completely static during the electron spin time evolu- 
tion. This is justified because the correlation time of 
the fluctuations in the nuclear-spin system due to dipole- 
dipole and hyperfine-mediated interaction between the 
nuclear spins, which is predicted to be > 10 — 100 (is 
IE, EI El ES EE El, III, HI, is much larger than the 
timescale for electron spin dynamics considered here (up 
to 1/is). 

In the experiment, the electron spin state is detected 
in a regime where electron transport through the dou- 
ble quantum dot occurs via transitions from spin states 
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with one electron in each dot (denoted as (1,1)) to the 
singlet state \S(0, 2)) with two electrons in the right dot. 
These transitions, governed via the tunnel coupling t c by 
the tunneling Hamiltonian H tc = t c \S(l, 1)) (5(0,2)1 + 
H.c, are only possible for anti-parallel spins, be- 
cause (TT|fltJS(0,2)) = (U| H u |S(0, 2)) = 0, while 
<4TI i?t |S(0,2)> , <n| i?t o JS(0, 2)) 0. Therefore, the 
states with even spin parity (parallel spins) block trans- 
port, while the states with odd spin parity (antiparallcl 
spins) allow for transport. If the system is initialized to 
an even spin-parity state, the oscillating transverse mag- 
netic field (if on resonance) rotates one (or both) of the 
two spins and thus lifts the blockade [l8j]. Initializing to 
|t) in both dots (the case with ||) gives the same result), 
we calculate the probability for an odd spin parity P dd 
under time evolution for each of the two spins governed 
by the Hamiltonian in Eq.{T]). 

Introducing the detuning from resonance 8 U = t z + 
h z — lo, the probability to find spin up for a single value 
of h z in the rotating wave approximation (which is valid 
for (b/e z ) 2 <C 1) is given by 
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b 2 + A8 2 b 2 +45l 

(2) 

Assuming that u> = /in + e z , i.e., 8^ = h z — /in, we find 
when av erag ing over the Gaussian distribution of h z val- 
ues (see [24[) 
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> max(i,l/&,6/2cr 2 ), with C = \ - 
' erfc I — S — ] . We can now calculate the 
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probability of finding an odd spin-parity state taking 
u> = ho + e z for both dots and drawing the value of h z 
independently from a distribution with width a in each 
dot: 



PoddW = P tL (t)(l - P tR (t)) + (1 - P hL (t))Pt,R(t) 



/(*) g(t) 



fit) 
g(t) 
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This result is valid for times t > max(l/cr, 1/6, b/2a 2 ) ~ 
20ns for a 1.4 mT nuclear field (see below) and b < 2a 
(accessible experimental regime). The 1 /i-term oscillates 
with the double Rabi frequency which is the result of both 
spins being rotated simultaneously (see also [HI]). This 
term only becomes important for b > a, because in that 
case for both spins most of the nuclear-spin distribution 
is within the Lorcntzian lincshape of the Rabi resonance. 



The l/v^-term oscillates with the Rabi frequency and 
orig inates from only one of the two spins being rotated 
[181 ] . This term is important when b < a , i.e., when only 
a small fraction of the nuclear-spin distribution is within 
the lineshape of the Rabi resonance. 

We also give the expression for P dd(t) f° r the case 
where only one of the two spins is on resonance (e z + ho — 
uj = 0), while the other is far off-resonance (\e z +ho~- u>\ 3> 
a). In this case the spin in one dot always remains up 
while the spin in the other dot rotates. This leads to 
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(7) 



with the same range of validity as in Eq. (f4|) . We see that 
the 1/t-term, which oscillates with frequency &, is not 
present in this case. 

The expressions for P dd(t) (Eqs. (f4]) and 0) reveal 
two interesting features: the power-law decay and a uni- 
versal phase shift of ir/4 (see Eq. ([5])) in the oscillations 
which is independent of all parameters. These features 
can both only appear if the nuclear field h z is static dur- 
ing a time much longer than the Rabi period. This is 
■ crucial because only then the driven spin coherence for 
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FIG. 1: (Color online) Rabi oscillations for four different driv- 
ing fields Sac (B z =55 mT, g=0.355 and a = gp, B {lA mT)). 
The gray circles represent the experimentally measured dot 
current (averaged over 15 s for each value of t), which re- 
flects the probability to find an odd spin-parity state after the 
RF burst that generates B ac . The dotted, solid and dashed 
lines represent the best fit to the data of an exponentially de- 
caying cosine function and the derived analytical expressions 
for Poddit) and P^ dd {t) (Eqs. © and ©) respectively. For 
clarity, the dashed line is shown only for the top two pan- 
els. The fit was carried out for the range 60 to 900 ns and 
the displayed values for _B ac were obtained from the fit with 
Podd(t) (Eq. ©)• We fit the data with an exponentially de- 
caying cosine with a tunable phase shift that is zero at t = 0: 
aie~ t/a2 [cos{<f) - cos(27rf/a 3 + <j>)] + a 4 (l - e~ t/a2 ). The last 
term was added such that the saturation value is a fit pa- 
rameter as well. We note that the fit is best for <j> — 7r/4, as 
discussed in the text. 
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FIG. 2: a) Dot current after an RF burst of 950 ns as a 
function of B ac , approximately representing the steady-state 
value. The solid curve is the best fit with ai(^ — 2C 2 ): 
the steady state expression of Eq. Q with ai and a as 
fit parameters. We find, for the 95%-confidence interval, 
a — g^is(1.0 — 1.7 mT). b) Decay power obtained from the 
best fit of the data (partially shown in Fig. [1} with the ex- 
pression ai + a2 cos(27ri/a3 + 7r/4)/t d , where 01,2,3 and d are 
fit parameters. 



one fixed value of h z is fully preserved. Because differ- 
ent values of h z give different oscillation frequencies, the 
decay is due to averaging over the distribution in h z . 

The phase shift is closely related to the power-law de- 
cay because it also finds its origin in the off-resonant 
contributions. These contributions have a higher Rabi 
frequency and shift the average oscillation in phase. This 
universal phase shift therefore also characterizes the spin 
decay, together with the power law. Interestingly, the 
specific shape of the distribution in h z (as long as it is 
peaked around the resonance) is not crucial for the ap- 
pearance of both the power-law decay and the phase shift 
[2J]. The values of the decay power and the phase shift 
are determined by the dependence of the oscillation fre- 
quency on h z (in this case ^Jb 2 + 4<5J). 

A power-law decay has previously been found theo- 
retically in 13, 12, 2H, 2(| and both a power-law decay 
(1/i 3 / 2 ) and a universal phase shift also appear in double 
dot correlation functions 13|, 21 1 . In [l~7| a singlet-triplet 
correlation function was measured, but the amplitude of 
the oscillations was too small for the phase shift and the 
power-law decay to be determined. Here, we consider 
driven Rabi oscillations of a single electron spin with a 
power-law decay of 1/ \fi that is already valid after a 
short time 1/er ~ 20 ns. Therefore, the amplitude of the 
driven spin oscillations is still high when the power-law 
behavior sets in, even for small driving fields (b < la) 
which are experimentally easier to achieve. The power- 
law decay and the phase shift thus should be observable 
in the experiment. 

We now discuss the observation of the power-law decay 
in the experimental data of which a selection is shown in 
Fig. [T] The data are obtained with the same device 
and under the same experimental conditions as in 
A fit is carried out to the observed oscillations for four 
different driving fields B ac (Fig. [T|), with three different 
fit functions: the theoretical expressions (Eqs. (HJ and 



([7]) with b and a constant scaling factor as fit parameters) 
and an exponentially decaying cosine. The width of the 
nuclear distribution a = gfj,g(lA mT) is obtained from 
a fit of the steady state value \ — 2C 2 of P dd{t) to a 
dataset obtained at t = 950 ns (Fig. [5h.). 

For the range B ac > 1.9 mT, we find good agreement 
with the model that predicts a power-lay decay of 1/y/i 
(Eq. ([4]); ho equal for both dots), while the fit with an 
exponentially decaying cosine is poor (blue lines in Fig. 
HJ. The power of the decay is independently verified by 
means of a fit to the data with ai+a 2 cos{2irt / a^+iv / A) / t d 
where, besides 01,2.3, the power d of the time t is a fit 
parameter as well. We find values of d ~ 0.6 (Fig. [2Jd), 
close to the predicted 1/Vt-dependence. 

We see much better correspondence of the data with 
Eq. (j4j than with Eq. 0, from which we can conclude 
that the mean of the Gaussian distribution ho is compa- 
rable for both dots (in equilibrium, we expect ho ~ in 
both dots). There might however still be a small differ- 
ence in ho between the two dots, which we cannot de- 
termine quantitatively because the two models describe 
only two limiting cases. If present, such a difference in 
ho could help explain the small deviation between data 
and model at the first oscillation for B ac = 2.5 mT. It 
could originate from asymmetric feedback of the electron 
spins on the respective nuclear spin baths, e.g. due to 
unequal dot sizes, leading to different hyperfine coupling 
constants. 

Another observation is that for small driving fields, 
B ac < 1.9 mT, we see that the damping is faster than 
predicted. Possible explanations for this effect are correc- 
tions due to electron-nuclear flip-flops (transverse terms 
in the hyperfine Hamiltonian) or electric field fluctua- 
tions. Electron-nuclear flip-flops may become relevant on 
a timescale ~ e z /cr 2 ~ 1 fj,s in this experiment. Electric 
field fluctuations can couple to spin states via the spin- 
orbit interaction or a finite electric-field dependent 
exchange coupling. 

We continue the discussion with the experimental ob- 
servation of the second theoretically predicted prominent 
feature of the Rabi oscillations, i.e., a phase shift of 7r/4 in 
the oscillations, which is independent of all parameters. 
The value of <p can be extracted most accurately from the 
oscillations measured for a wide range and small steps of 
B ac . like the data shown in Fig. [3^,. That is because the 
Rabi period T Rabl = 2tt/ g^ B (^B ac ) = 2n / g^ B {\KI s ) 
contains only one unknown parameter K (current to os- 
cillating field amplitude B ac conversion factor, in units 
of T/ A) which is independent of the current through the 
wire I s that generates B ac [l8[ . The presence of a phase 
shift is visible in Fig. [3k,, where the green and blue lines 
correspond respectively to the maxima of a cosine with 
and without a phase shift of 7r/4. The green lines match 
very well the yellow bands representing high data values. 
In contrast, the blue lines are located on the right side 
of the yellow bands for small burst times and more and 
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FIG. 3: a) The dot current (represented in colorscale) is 
displayed over a wide range of B a c (the sweep axis) and burst 
durations. The green and blue lines correspond respectively 
to the maxima of a cosine with and without a phase shift of 
7r/4. The current-to-field conversion factor K is fitted for both 
cases separately (A"=0.568 niT/mA and 7^=0.60 mT/mA for 
respectively with and without phase shift; the fit range is 
t = 60 - 500 ns and I B = 3.6 - 6.3 mA). b) Phase shift 
for a wide range of _B ac , displayed as a function of stripline 
current I s . Values obtained from a fit of each trace of the 
data in a) (varying burst time, constant -Bac) to a damped 
cosine a\ — ai cos(^KI B g/j,Bt + as^/y/t, where ai,2,3 are fit 
parameters and K — 0.568 mT/mA. I s is a known value in 
the experiment, extracted from the applied RF power. The 
gray dashed lines represent the 95% confidence-interval. 



more on the left side of the bands for increasing burst 
times. Thus, a cosine without a phase shift does not 
match with the observed Rabi oscillations. 

In order to determine <f> quantitatively, we perform 
a single two-dimensional fit of the complete dataset in 
Fig. with P dd{t) (Eq. fl2J)), excluding the 1/t-term 
(see Hi). The fit range is t = 100 - 900 ns, such 
that the contribution from the 1/t-term of Eq. ([4]) can 
be neglected. For the 95% confidence interval we find 
<j> = (0.23 ± 0.01)7r, close to the theoretical value. The 
relation between (j> and B ac is visible in Fig. where we 
find no significant dependence of <f> as a function of -B ac , 
although the accuracy decreases for smaller B ac (values 
obtained from fits to single traces, see caption). We have 
not compensated for the effects of the finite rise time (<2 
ns) of the bursts, which leads to a small negative phase 
shift, on top of the expected positive 7r/4 shift. 

To conclude, we have experimentally observed a power- 
law decay and universal phase shift of driven single elec- 
tron spin oscillations. These features are theoretically 
understood by taking into account the coupling of the 
spin to the nuclear spin bath, which is static on the 
timescale of the electron spin evolution time. Further- 
more, the slow power-law decay allows spin manipula- 
tion with relatively small driving fields. This improved 
understanding of the coherence of a driven single elec- 
tron spin is important for future experiments using the 
electron spin as a qubit. For future investigation, it re- 



mains interesting to obtain more information about the 
non-static contributions of the nuclear bath or other pos- 
sible decoherence mechanisms. For that, it is required to 
measure the driven oscillations at larger external fields, 
with larger driving powers and longer evolution times 
than accessible in this work. 

We thank T. Meunier, R. Hanson, Y.V. Nazarov and 
I.T. Vink for discussions; R. Schouten, A. van der En- 
den, R. Roeleveld and W. den Braver for technical as- 
sistance. We acknowledge financial support from the 
Dutch Organization for Fundamental Research on Mat- 
ter (FOM), the Netherlands Organization for Scientific 
Research (NWO), JST ICORP, NCCR Nanoscience, and 
the Swiss NSF. 



[1 
[2 

[3] 

[4] 

[5 

[6 

[~\ 
[« 
[9 
[10 

[11 

[12 

[13 

[14 
[15 
[16 
[17 
[18 
[19 

[20 

['21 

[22 
[23 
[24 

[25 

[26 

[27 



Rev. 



Rev. 



F. Bloch, Phys. Rev. 70, 460 (1946). 

D. P. DiVincenzo and D. Loss, Phys. Rev. B 71, 035318 
(2005). 

J. Taylor and M. Lukin, Quantum Information Process- 
ing 5, 503 (2006). 

G. Ithier et al, et al., Physical Review B 72, 134519 
(2005). 

A.V. Khaetskii and Y.V. Nazarov, Phys. Rev. B 61, 
12639 (2000). 

V. N. Golovach, A. Khaetskii, and D. Loss, Phys. 
Lett. 93, 016601 (2004). 

J. M. Elzerman et al, Nature 430, 431 (2004). 
M. Kroutvar et al, Nature 432, 81 (2004). 
S. Amasha et al, arXiv:cond-mat/0607110 (2006). 
A. V. Khaetskii, D. Loss, and L. Glazman, Phys. 
Lett. 88, 186802 (2002). 

LA. Merkulov, A.L. Efros, M. Rosen, Phys. Rev. B 65, 
205309 (2002). 

W. A. Coish and D. Loss, Phys. Rev. B 70, 195340 

(2004) . 

W. A. Coish and D. Loss, Phys. Rev. B 72, 125337 

(2005) . 

A. C. Johnson et al, Nature 435, 925 (2005). 

F. H. L. Koppens et al, Science 309, 1346 (2005). 

J. R. Petta et al, Science 309, 2180 (2005). 

E. A. Laird et al, Phys. Rev. Lett. 97, 056801 (2006). 

F. H. L. Koppens et al, Nature 442, 766 (2006). 

W. Yao, R.-B. Liu, and L. J. Sham, Phys. Rev. B 74, 
195301 (2006). 

R. de Sousa and S. Das Sarma, Phys. Rev. B 67, 33301 
(2003). 

D. Klauser, W. A. Coish, and D. Loss, Phys. Rev. B 73, 
205302 (2006). 

C. Deng and X. Hu, Phys. Rev. B 73, 241303(R) (2006). 
C. Deng and X. Hu, Phys. Rev. B 74, 129902(E) (2006). 
See EPAPS Document No. [number will be inserted by 
publisher] for fit procedures and asymptotic expansion. 
V. V. Dobrovitski, H. A. De Raedt, M. I. Katsnelson, 
and B. N. Harmon, Phys. Rev. Lett. 90, 210401 (2003). 
R.-S. Huang, V. Dobrovitski, and B. Harmon, 
arXiv:cond-mat/0504449 (2005). 

M. Borhani, V. N. Golovach, and D. Loss, Phys. Rev. B 
73, 155311 (2006). 



5 



SUPPLEMENTARY MATERIAL FOR 
"UNIVERSAL PHASE SHIFT AND 
NON-EXPONENTIAL DECAY OF DRIVEN 
SINGLE-SPIN OSCILLATIONS" 

Fit Procedures 

Here, we describe the exact procedure of the two- 
dimensional fit from which the phase shift <j> = (0.23 ± 
0.01) was obtained. The fit function is a simplification 
of -Podd(i) (Eq. (4)). The first simplification is the ex- 
clusion of the 1 / t-term because its contribution is negli- 
gible within the fit range of t = 100 — 900 ns. Second, 
both expressions for C and 0.5 — 2C 2 are approximated 
as being linear in b, which is justified for the regime we 
consider (I s = 3.6 — 6.3 mA), as can be seen in Fig. 
[TJ With these simplifications, we obtain the expression 
2/modoi = a 1 +I s a 2 + {a 3 +I s a4)cos(^KI s gfi B t/h+(/)Tr)/\/i 1 
where 01,2,3,4, K and (f> are fit parameters. The Rabi fre- 
quency WR a bi is given by \KI s g\iB /h, with I s the current 
through the stripline which is known in the experiment. 
The constant factor K = B ac /I s is not known in the ex- 
periment but can be obtained from the fit. The behavior 
of J/modei — J/data around the optimal values for the fit 
parameters is seen in Fig. [2] 

As a cross-check, we carried out a fit for each 
trace of the data in Fig. 3a of the main text (vary- 
ing burst time, constant £? ac ) with a damped cosine 
di — a-i cos(^KI s gfiBt/h + <jm)J^/t, where a\ t 2 and <f> are 
fit parameters, and K is kept at a constant value. This 
fit is carried out for a wide range of values for K and 
I s . The best- fit values for <j> obtained for different I s arc 
averaged and plotted as a function of K in Fig. 03 (blue 
curve), together with the spread in 4> (gray dotted lines) 
and the fit quality (green curve). The figure shows that 
the best fit is obtained for </> = 0.237T, the same value we 
found from the two-dimensional fit. 
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FIG. 2: Phase shift and fit quality for a range of current-to- 
field conversion factors K. Values obtained from a fit with a 
damped cosine to the single traces (constant I s ) of the data set 
shown in Fig. 3a of the main text, and subsequently averaged 
over all traces for the range t — 100 — 900 ns and I 3 =3.6- 
6.3 mA. The fit-quality R 2 is a measure of the correlation 
between the observed values t/i and the predicted values 

* 2 =Er =1 i5f#>withy=i£^ 
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FIG. 3: Root mean square difference between 

the measured current j/data and the model y mo del : 

fy~^ (Vtff - y"° a de ') 2 ^ /N t Ni,, for a wide range of 

K and <j>. We sum over the range t = 100 — 900 ns and 
I s =3.6-6.3 mA. 



Asymptotic expansion 

Here we give steps and additional justification leading 
to the asymptotic expansion given in Eq. (3) of the main 
text. We consider averaging Eq. (2) from the main text 
over a quasicontinuous distribution of h z values for the 
case where 6 U = h z — ho (with the replacement 8 U — > x) : 



FIG. 1: C and 0.5 — 2C 2 as a function of I s , with a = <?^b(1.4 
mT) and b = g^isKIs- K = 0.56 mT/mA for the curves 
shown here, but the curves remain linear as well for if=0.5- 
0.6 mT/mA. 



P T (f) = / dxD{x)P hx {t). (1) 
J — 00 

As a consequence of the central-limit theorem, for a large 
number of nuclear spins in a random unpolarizcd state, 
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the distribution function D(x) is well-approximated by a 
Gaussian with standard deviation a centered at x = 0: 



D(x) 



1 



2?ra 



(2) 



Inserting Eq. ([2|) into Eq. {1} gives the sum of a time- 
independent part, which can be evaluated exactly, and a 
time-dependent interference term I(t): 



(3) 



Here, C is given following Eq. (3) of the main text. With 
the change of variables u = (V& 2 + 4a; 2 —6) /2a, and 
using the fact that the integrand is an even function of 
x, the interference term becomes I(t) = Rel(t), where 
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When at ^> 1, the time dependence of I(t) is controlled 
by the region u < 1/at. The integrand simplifies con- 
siderably for au/b <C 1, which coincides with at ^> 1 
for 



z 1 b 1 
it < — <C i > -. 
at a a 



Equivalcntly, for 



t > max 



1 1 

Va 



(5) 



(G) 



(4) we expand the integrand for u < min(l, &/ct): 



I(t) 



Jbt/2 



&ira 



(III 



exp | 



-Xu + 0(u 2 )) 



1 + 



au\ 

~b) 



r 



b 

2^ 



(7) 



Neglecting corrections of order u 2 in the exponential and 
order au/b in the integrand prefactor, the remaining in- 
tegral can be evaluated easily, giving 



1(f) 



cos [bt/2 + arctan(i/r)/2] 



t > max (1/6, 1 /a) . 



1/4 



-, T 



b/2a 2 , (8) 
(9) 



Eq. ([S]) is valid for the time scale indicated for an ar- 
bitrary ratio of b/2a. Due to the exponential cutoff at 
u < 2a /b in Eq. 0, Eq. © is actually valid for all 
times in the limit b/2a~^> 1. Expanding Eq. ([8]) to lead- 
ing order for t/r 1 gives the result in Eq. (3) of the 
main text. Higher-order contributions to the long-time 
expansion of Eq. ([8]) and contributions due to correc- 
tions of order au/b in Eq. ([7]) both lead to more rapidly 
decaying behavior of order ~ 1/i 3 / 2 . The reason for the 
different phase shift here (7r/4) relative to that found in 
Ref. [l3| (37r/4) is that here the fluctuations are longi- 



tudinal, while in Ref. [13| the fluctuations are along the 
transverse direction. This leads to a different integrand 
in Eq. (S4) and thus to a different value for the phase 
shift and decay power. 

Since the long-time behavior of I(t) is dominated by 
the form of the integrand near x = 0, the same result 
can be found after replacing the Gaussian distribution 
function by any other distribution function D(x) which 
is analytic and has a single peak at x = 0. Specifically, 



D{x) 
1 



D(0) exp{- 
d 2 IylD(x) 



2a 2 



0(x 3 )}, 



dx 2 



(10) 
(11) 



x=0 



Thus, the universal form of the long-time power-law de- 
cay and phase shift are relatively insensitive to the spe- 
cific shape of the distribution function. 



